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Abstract 

We consider models of accelerating Universe elaborated for Finsler 
like gravity theories constructed on tangent bundles to Lorentz man- 
ifolds. In the osculating approximation, certain locally anisotropic 
configurations are similar to those for f(R) gravity. This allows us 
to generalize a proposal (by Nojiri, Odintsov and Saez-Gomez, AIP 
Conf. Proc. 1458 (2011) 207-221; arXiv: 1108.0767) in order to re- 
construct and compare two classes of Einstein-Finsler gravity, EFG, 
and f(R) gravity theories using modern cosmological data and realistic 
physical scenarios. We conclude that EFG provides inflation, acceler- 
ation and little rip evolution scenarios with realistic alternatives to 
standard AC DM cosmology. The approach is based on a proof that 
there is a general decoupling property of gravitational field equations 
in EFG and modified theories which allows us to generate off-diagonal 
cosmological solutions. 

Keywords: modified theories of gravity, Finsler geometry, accel- 
erating Universe, dark energy and dark matter. 

PACS: 04.50.Kd, 04.20. Cv, 98.80.Jk, 95.36. +x, 98.80.Cq 



*pstavrin@math. uoa.gr 

t sergiu.vacaru@uaic.ro, http://www.scribd.com/people/view/1455460-sergiu 



1 



1 Introduction and Preliminaries 



Modern cosmology is based on observational data for two accelerating 
periods in evolution of Universe, the early-time inflation phase and the late- 
time acceleration with dark energy and dark matter effects. Such acceler- 
ating epochs are characterized by a number of differences with large and, 
respectively, small values of curvature; quite exotic states of matter at the 
"beginning" and "end"; alternative classes of solutions of the gravitational 
and matter fields equations, with singularities and possible anisotropies etc. 
One of the main exploited ideas is that the Universe is under cyclic evolution 
with oscillating equations of matter [TJ [2J O 0] . It is considered that it may 
exist an unified theory which describes in different limits both the inflation 
and dark (energy/matter) periods. 

Modified gravity theories were elaborated as unifications and/or gener- 
alizations to Einstein gravity when cyclic cosmology with inflation and dark 
energy may be realized (see reviews [U EJ [7]). The approaches with mod- 
ified Lagrange density R — > f(R,T, ...), where R is a scalar curvature and 
a functional /(...) is determined by traces of certain stress-energy tensors, 
additional torsion fields etc, are intensively developed in modern literature. 
Nevertheless, there are other directions which at first site are alternative to 
/(...) gravity and present interest in modern classical and quantum gravity 
and cosmology. In this work, there are studied the so-called Finsler gravity 
and cosmology (for review of results and critical remarks, see [8| l9l fTUl [TT] ) 
and investigate possible connections to f(R) theories. Perhaps, the first 
model with locally anisotropic inflation in a Finsler like manner was pro- 
posed in [12]. Several Finsler cosmology and gravity models where devel- 
oped in Refs. [ia[I3[l5l[lS[l7l[ll[ll[20l[2ll[22] following geometric and 
physical ideas related to quantum gravity and modified dispersion relations, 
broken local Lorentz symmetry, nonlinear symmetries etc. 

Our purpose is to prove that a class of Einstein-Finsler gravity (EFG) 
theories can be considered as natural candidates for which cyclic cosmology 
with inflation and dark energy/matter epoches can be realizedQ We shall 
formulate a reconstruction procedure when alternatively to f{R) gravity 
|23j we can extract information on certain locally anisotropic (Finsler) like 
gravitational models. It should be noted that the general relativity (GR) 
theory and modifications can be re-written in so-called Finsler like variables 

1 An EFG model is constructed on a manifold V, or its tangent bundle TV, similarly to 
the GR theory when, roughly speaking, the Levi-Civita connection V is substituted by a 
Finsler like metric compatible connection D, both defined by the same metric structure 
g. The second linear connection D is with nontrivial torsion. 
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which allows us to decouple and integrate the gravitational and matter field 
equations in very general forms. This is possible for nonholonomic/non- 
integrable 2 + 2 splitting of a (pseudo) Riemannian manifolds enabled with 
a formal fibred structure. We can mimic a (pseudo) Finsler geometry on 
a standard Einstein manifold if we adapt the constructions with respect 
to corresponding nonholonomic (non-integrable) distributions. Such con- 
structions are rather formal but allow us to formulate a general geometric 
method of constructing exact solutions both in GR and f(R, T, ...) modifica- 
tions, when the generic off-diagonal metrics and various types of connection 
and frame variables depend on all coordinates via generating and integration 
functions and parameters. 

Let us summarize some key ideas on modifications/ generalizations of the 
Einstein gravity on (co) tangent bundles. In such theories, the fundamental 
geometric objects (metrics, frames and connections) depend on velocity/ 
momentum type variables which can be interpreted as fiber like coordinates. 
The geometric constructions are derived from a nonlinear quadratic element 
ds 2 = F 2 (x,y), where x = (x l ) and y = (y a ) are coordinates on a tangent 
a bundle TV to a Lorentz manifold V. Such base manifolds are necessary 
if we wont to get in a limit the standard GR theoryU The value F is called 
the fundamental/generating Finsler function. Usually, certain homogeneity 
conditions on y are imposed and the nondegenerated Hessian 

~ 1 d%F m 
9ab 2 dy a dy b 1 j 

is considered as the fiber metric. Gravity theories with anisotropics on TV 
can not be determined only by F, or g a ^. This is very different from GR 
and a class of f(R,T, ...) theories (on pseudo-Riemannian spaces) when the 
geometric and physical models are completely defined by data the metric 
structure. To construct a Finsler geometry/gravity theory we need a triple 
(F : g, N,D) of fundamental geometric objects: the total metric structure, 
g, the nonlinear connection (N-connection) structure, N, and the distin- 
guished connection structure, D, which is adapted to N. There are necessary 
additional assumptions how such objects are defined by F, which types of 
compatible or noncompatible linear connections are involved, how the cor- 
responding curvature, torsion and nonmetricity fields must be computed. 
For realistic physical models, the experimental/observational effects can be 



2 The well known pseudo-Riemannian geometry consists a particular case with 
quadratic element ds 2 — gij(x)dx' 'dx\ when y 1 ~ dx l . 
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analyzed via a so-called osculating approximation (see examples in |13|,I14|) 

1 d 2 F 

9i j =g ij (x,y(x)) = 2d^b( X ' y ( X V- (2) 

For a review of Finsler geometry for physicists, see Refs. [H [11] and refer- 
ences therein. 

The paper is structured as follows. In section [21 we summarize the 
necessary results on Einstein-Finsler and f{R) modified gravity and show 
how such theories can modelled by nonholonomic distributions and/or off- 
diagonal metrics in GR. Section [3] is devoted to models of conformal cyclic 
universes in EFG. We consider FLRW metrics subjected to nonholonomic 
constraints and analyze consequences of locally anisotropic models. A pro- 
cedure of reconstructing EFG theories is provided. There are studied several 
models with effective anisotropic fluids. Scenarios of ekpyrotic and little rip 
cosmology governed by fundamental Finsler functions in osculating approx- 
imation are studied in section SI The anholonomic deformation method is 
applied for generating cosmological solutions in EFG. Conclusions are drawn 
in section [5j 

2 Canonical EFG & f(R) Modifications 

We consider a four dimensional (4— d) Lorentz manifold V in GR mod- 
elled as a pseudo-Riemannian space enabled with a metric ft g = {gij{x 1 )} 
of local signature (+, +, +, — ). Physically motivated Finsler generalizations 
to metrics and other geometric objects depending anisotropically of veloc- 
ity/momentum type variables y a can be constructed on tangent bundle TV. 
A (pseudo) Finsler geometry is characterized by its fundamental (equiva- 
lently, generating) functions when certain homogeneity conditions are im- 
posed, F(x l 1 f3y : ') = f3F(x l ,y 3 ), for any (3 > 0, and det \g a b\ 7^ 0, see Hessian 
(PQ) which defines a so-called vertical metric. In standard approaches, the 
matrix g a ^ is considered positively definite, but this condition has to be 
dropped in (pseudo) Finsler geometry (hereafter, we shall omit the term 
"pseudo"). 

2.1 The "triple" of fundamental geometric objects 

A canonical model of Finsler-Cartan geometry is completely determined 
by F, and gij, up to necessary classes of frame/coordinate transform e a = 
e a a (x,y)e a , following such assumptions on the triple (F : N,g, D) of funda- 
mental geometric objects (see details in [8| fTTj): 
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There is a canonical nonlinear connection (N-connection) structure, 



N : TTV = hTV vTV, (3) 

which can be introduced as a nonholonomic (non-integrable/ anholonomic) 
distribution with horizontal (h) and vertical (v) splitting]! We get an inte- 
grable/holonomic frame configuration if W^g = 0. Under frame transforms, 
the coefficient formulas transform into equivalent ones for arbitrary sets of 
coefficients N = {N", = e^ejN?}. We can adapt the geometric construc- 
tions via " N-elongated" (co) frame structures, 

e u = ( e ,,e a ), e i = A _iVf(u)A ^ 6a = _^_ ; ( 4 ) 

= (e', e a ), e* = dx i and e a = dy a + N?{u)dx\ (5) 
In general, such frames are nonholonomic because 

[e a , e^] = e a e/3 - e^e^ = Wj /3 w 7 (6) 

i i r)N a ON? i.dN a i.f)?J a 

with W\ a = d a N> and W £ = 0% = + - N*^. 

A canonical metric structure on TV, can be introduced using data 
(9ij,e a ) , with e a > = e a a ,e a and g a /p> = e a a , eT^ g af3 , 

g = hg vg = gij (x, y) e* <g> e j + h ab (x, y) e a ® e b . (7) 

The third fundamental geometric object in a Finsler geometry is the 
distinguished connections (d-connection) D = {T^ 7 } =(hD,vD) which by 
definition is adapted to the N-connection structure, i.e. preserves the non- 
holonomic h-v-splitting ([3]). We can model physically viable "almost" stan- 
dard models (see discussions and critical remarks in [HE]) for d-connections 
which are compatible with the metric structure, Dg = and completely de- 
fined by data (g, N). Let us consider how such a canonical d-connection D 



3 To define the canonical N-connection we follow a geometric/variational principle for 

l 

an effective regular Lagrangian L — F 2 and action S(t) = J L{x(r),y{T))dr, for y k (r) — 

o 

dx k (r)/d,T. The Euler-Lagrange equations J^t^jt — J^r = are equivalent to the "non- 
linear geodesic" (equivalently, semi-spray) equations d d *. i +2G h (x,y) = 0. The canonical 
coefficients N = {N?} are computed Nf := ^1^1, G k = \~g h] (V - 

In our works, we put "tilde" on symbols (for instance, e a ) if the constructions are 
performed with respect to bases determined by N°" . "Tilde" will be omitted for general 
N-adapted h- ^-splitting when variables mix each to others. 



5 



can be constructedo By definition, such a connection is with zero h- and 
v -torsions, T l - k = and T a bc = 0. There are also nontrivial torsion coeffi- 
cients. In Finsler geometry, it is also possible to introduce in standard form 
the Levi-Civita connection V = {rjg } (as a unique one which is metric 
compatible and zero torsion) but it does not preserve under parallelism the 
N-connection splitting (|3|). There is a canonical distortion relation 

D =V + Z, (8) 

where D, V and the distortion tensor Z are uniquely defined by the same 
metric structure g. 

Finally, we note that there is an important argument to work with D 
with respect to N-adapted frames and ([5]) even in GR. This is possible for 
a conventional nonholonomic 2+2 spitting. The priority is that D allows to 
decouple and integrate the gravitational field equations in very general forms 
[201 IHl E2] (see section 0|) . Such solutions define standard Lorentz /Einstein 
manifolds if we constrain at the end the generating/integration functions in 
such a form that Z = which mean zero torsion. For such configurations, 
we get D| T=0 =V in N-adapted form, see (JHJ). 

2.2 The EFG field equations and f(R) gravity 

A Einstein-Finsler gravity theory (EFG) is constructed for a Finsler like 
d-connection D following standard geometric and/or variation rules as in 
GR. In this work, we shall use only the canonical d-connection. Such models 
can be constructed, for instance, on a Lorentz manifold V (considering an 
2+2 splitting; in this case, we introduce Finsler like variables in GR, when 
T = is considered as nonholonomic constraint for D — > V, or on its tangent 

4 Using N-adapted differential forms and the d-connection 1-form is r a s = r a p 1 e"' , 
we compute the torsion, T a , and curvature 2-forms, 1Z a n, T a := De Q = de a + 
A e' 3 , n a p := Dr a s = dC a p - A T a 7 = R^e 7 A e 6 . The torsion coefficients 

Ta f rpi rjii rya rjia rjia \ T i T * T 1 ^ T"^ rjia 

0-/ — I 1 jk; 1 jai 1 ji; 1 bi; 1 be I arc 1 jk — ^ jk ~~ u kj i 1 ja — ^ 1 aj — jai 1 ji — 

T a u = ^0r~ La bi, T a bc = C a bc - C a cb ; for N-adapted coefficients TL% S , see [8] 03]. 
In our former works the symbol D was used for the canonical d-connection. We shall 
omit "hats" if that will not result in ambiguities. The N-adapted coefficients of D can be 
computed in the form r 7 Q/3 = (L) k , L£ fe , C) c , C b a c ) , 

Ljk — 2^ ( efc fJ r + e i9kr — e r gjk) , C% c = -zh ad (e c h b d + e c h c d — edh bc ) , 
Ltk = e 6 (7Vfe) + i/i ac (e k h bc - h dc e b N% - h d b e c N^j , C) c = -g tk e c gjk- 
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bundle TV. In the second case, we can relate the constructions to certain 
models of Finsler-Cartan gravity (via nonholonomic deformations, we can 
transform D into the Cartan connection for Finsler space 

The Ricci tensor, Ric = {R ^ := R a ^ 7Q ,}, of D =(hD,vD), splits into 
certain h- and u-components (Rij , R a i, Ri a , R a b). The scalar curvature is 

F R := R^ = gVRij + h ab R ab = R + R. (9) 

The gravitational field equations can be postulated in geometric form (and/or 
derived via N-adapted variational calculus), 

R^-gS/w F R = T P s- (io) 

Such equations transform into the Einstein's ones in GR if T ps —> Tp§ and 
D — > V. So, we can consider necessary N-adapted variations of actions for 
scalar, electromagnetic, spinor etc fields in order to derive certain the matter 
source in (|lUp . The solutions of the field equations in EFG are with nontrivial 
torsion induced by (g, N). The Levi-Civita (zero torsion) conditions with 
respect to N-adapted frames are 

^• = e a (ivj),q 6 = o, n a Jt = o. (ii) 

If such constraints are imposed at the end after we have constructed certain 
classes of solutions of (|10p on a 4-d manifold and for physically motivated 
sources, we generate (in general, off-diagonal) solutions in GR. 

It is not clear what types of sources should be considered in EFG 
models on TV. We note that, in general, such tensors are not symmetric 
because R p$ is not symmetric. It reflects the nonholonomic and locally 
anisotropic character of Finsler gravity theories. For certain toy models, 
we can approximate to a cosmological constant with possible locally 

anisotropic polarizations depending on (x*,?/ a ) . The fundamental Finsler 
function F(x, y) is encoded into geometric objects of (flOl) and solutions of 
such PDE. In general, we can not "see" explicit dependencies on F(x,y) 
because of principles of generalized covariance and relativity: arbitrary 
frame/coordinate transforms on V and TV mix the variables. Fixing a 
system of local coordinates, we can measure experimentally on /i-subspace 
only an osculating (pseudo) Riemannian metric ([2]) . Mathematically, we can 

5 We use such conventions on indices a = (i, a) : For a 2+2 splitting on V, in GR, we 
consider i, j, ... = 1,2 and a, b, ... — 3, 4. On TV, ... — 1, 2, 3, 4 and a, b, ... = 5, 6, 7, 8. 
Here we note that formal integrations of gravitational field equations in general forms are 
possible for splitting of type 2 + 2 + 2 + 2 + ... and/or 3 + 2 + ... 
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construct exact solutions of (fTUj) for 8-d metrics (|7|) but to verify possible 
physical implications in direct form is possible only for the /i-components. 

EFG theories on V, or TV, are explicit examples of modified GR. In 
a different manner, such theories of gravity are constructed for modified 
Lagrange densities when R — > f(R, T, ...), see reviews of results in are [SJ El 
[7j. For instance, the corresponding field equations with f(R) can be written 
as effective Einstein fields equations 

iV - \R9w = &r ef G ef T^ u , (12) 

for ^T^ = (dRfr^lf-RdRfjg^-ig^Va-V^V^dRf] and = 

Using observational/ experimental data, we can measure some geometric 
configurations determined by functional dependencies of type gij[F, /,...]. In 
general, it may be not clear if such a metric is a solution of (I10p or (I12D . i.e. 
we can not say exactly what kind of modifications of GR would result in dark 
energy/matter effects. If a fundamental Finsler function F is considered as 
a nonholonomic distribution on a pseudo-Riemannian space V, or TV, mod- 
ifications of type R -> f(R) result in (F; F R) -»• f(R, F) ~ f(R) ~ /( F R), 
see formulas (jlOp for the Finsler curvature scalar and its h- and ^-splitting. 
In general, such modifications can be performed for any model of Finsler 
spacetime geometry with scalar curvature F R. A physically realistic theory 
closed to GR and MG can be constructed in the simplest way for the scalar 
curvature R of the Cartan d-connection. For a prescribed N-connection 
structure, we can define via nonholonomic deformations ([8]) functional de- 
pendencies of type R(R) and, inversely, R(R). In general, EFG and /(...) 
are different modifications of GR described by different Lagrange densities 
and derived field equations. For certain conditions, we can transform a class 
of solutions of (jlOp into (|12p . and inversely, using frame/coordinate trans- 
forms. It is known a reconstruction technique [23j allows us to recover data 
for f(R) using observational data from modern cosmology. In this work, we 
formulate a similar approach for extracting data for EFG theories and spec- 
ulate on conditions when we can distinguish a Finsler configuration from a 
/-modification. 
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3 Conformal Cyclic Finsler Like Universes 



In this section, we show how to construct cycling universes in models 
with f(R,F) ~ f(R) when the theory is determined by actions of type 



S 



d 4x \f\9ij\ 



f(R) 

2k? 



+ m L 



and/or S 



\9ij\ 



R(R) 
~2kJ~ 



+ m L 



(13) 

where k is the gravitational constant, m L is the Lagrangian for density and 
the osculating approximation (j2j) is encoded in gij,R and (for simplicity, 
we shall consider models with functional dependence on metric for matter 
Lagrangians) m L. 

A scalar field x can be introduced via conformal transforms g a n —> 



-x(z) 



g a /3 with respect to N-elongated frames (jlj and 1^1 Denoting 
X R = e*{R + 3g«[(D iX )(D iX ) - l(e iX )(e jX )}}, 



we get a new action 



S 



d A xe- 2 * 



9ij\ 



2k- 



:f( x R,X)+ m L 



(14) 



Choosing some "proper" functionals P(x) and Q(x) when the equation 

d x P( X ) R + d x Q( X ) = (15) 

can be solved as x = x(-^)> we can express f(R) = P[x(i?)] R + Q[x(R)]- 
For "pure" Finsler modifications, the functionals (|15p are of type R = 
R(R) with determined in nonholonomic form via distortions ([8]), where R is 
the scalar curvature of V. In such cases, actions of type (fl"4"|) are parameter- 
ized in the form 

1 

when 



S 



I 



d xe 



-2* 



9ij\ 



773 R ( R >X) + 



(16) 



R(R) = P[ x (R)}R + Q[x(R)}- (17) 

We conclude that an EFG theory with osculating approximation is a 
variant of modified gravity with / — > R. Such theories can be studied by 
the similar methods if we work with respect to N-adapted frames. In former 
approaches to f(R,T) modifications of gravity, such "preferred" bases and 
geometric constructions were not considered. 



We do not write <j> for scalar fields as, for instance, in [23] but use x because in 
our works cj> is considered as a generating function for constructing off-diagonal solutions 
[20l[2T][22] (see also next section). 
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3.1 Reconstructing EFG theories 

There are three possibilities to introduce locally anisotropic scalar fields 
in a iZ-modified theory: 1) via conformal transforms as in f)14[) : 2) with 
certain compactification of extra dimensions of models on TV; 3) postulating 
scalar field Lagrangians with operators ej and the canonical d-connection 
Dj. Models of type l)-3) mutually transform from one into another one 
under frame/coordinate transforms. 

Let us introduce a system of local coordinates (x 1 ,^ 2 ,^ 3 = i,x 4 ), with 
time like coordinate t, on V, and osculating approximation of Finsler type 
9ij = 9ij(t,y a (t)) and Nf = Nf{t,y h {t)) with "non-tildes" with respect to 
arbitrary frames. On /i-subspace, we can consider a FLRW ansatz of type 
(jZJ) when 

ds 2 = -(e 3 ) 2 + a 2 (t + t) [(e 1 ) 2 + (e 2 ) 2 + (e 4 ) 2 ] + {unknown v-components} , 

(18) 

with 

S(t) =exp[H t + b(t)], (19) 

where the constant Hq > and t is the period. We use tilde on a in order 
to emphasize that this value contains certain Finsler information. Such 
a generic off-diagonal metric is considered to result in physical identical 
scenarios at t and t + t. For instance, we can take b(t) = °6cos(27ri/i) 
when the condition °b < Hq results in monotonic expansion. We emphasize 
that the h— part of quadratic element on TV is the same as in the 4— d 
flat cosmology. Nevertheless, such models are with osculating dependence 
on t,y a (t) and contain certain information on a possible Finsler spacetime 
structure with N-elongated partial derivatives of type (HJ). To show this 
we consider a locally anisotropic scalar-tensor model. The action is taken 
following the possibility 3) above, 



S 



(20) 



for some functionals uj(x) an d V(x) on X- If both values are determined by 
a single function s(x), we express 

u( X ) = -(2/K 2 )d 2 xx ,( X ) and V( X ) = «" 2 {3[%(x)] 2 + d lx<X)} 

for a cosmological solution with h— part of FLRW type, when x = t and the 
Hubble function H = 9^0 



7 We note that for cosmological models, the rule of N-elongation of partial derivatives 
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For instance, if we parameterize the cycling factor (|19p a(t) = exp[<j(x)] = 
exp[<j(i)], we recover a scalar-tensor model with spherical symmetry on hTV 
and reproduce the cyclic universe via formulas which are similar to that pre- 
sented in Sec. II of [23], 



u( X ) = (2°6/ K 2 )(2^- 1 ) 2 cos(2vr X f- 1 ), 

V(x) = k- 2 {3[H - (2vr °bT 1 ) sin(2vrx^ 1 )] 2 - (27T?" 1 ) 2 cos(2vrx^ 1 )}. 

It is not possible to say to much on underling Finsler structure using such 
solutions. We can only conclude that a model with h-scalar R may be with 
cyclic behaviour. In order to extract more information, we should perform 
a rigorous analysis how F modifies the GR. 

We apply the methods of reconstructing modified gravity theories [H 
El [Ml E3] reformulated for models of type (fT3"]) with R(R) and functionals 
P[x]i? + Q[x] (|17p and scaling factor parameterized in the form a(t) = 
= const. The solutions for P and Q can be found expanding such 
values in Fourier series and defining certain recurrent formulas for coefficients 
expressed in terms of Ho and °b. For simplicity, we provide the equation for 
P when the matter can be neglected, 



^(x)/^/\^\\=-\[d xx P(x)}/(y/^\Y 

The equation for w(x) is 

d x ^ = — VIPI / d X (d xx P)/(x/\P\f. (21) 




We can consider a functional P in (|15p for certain values of R as a generating 
function for solutions of (|21j) . Here we provide a well known example with di- 
vergences corresponding to a(to) = which can be identified with a moment 
of singularity for a big bang, or crunch effect. Fixing P = Pq [cos(Pix)] 4 , 
for some constants Po and Pi, we express the last equation as 

3 x w = w [cos(Pix)] 2 + 2Pi [sin(2Pix) - tan(2P lX )] , 



0J for A(x(t),y(x(t))) is 

8 A 

d t A^ et A = d t A-N?^—, 
oy a 

which points to possible contributions from the fibred structure. Considering a scal- 
ing factor a(t) in a nontrivial nonholonomic (Finsler) background, we introduce H(t) := 
[e t o(t)] /a(t). 
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for an integration constant wq. The term with tan(2Pix) results in diver- 
gences of as it is well known for cyclic scenarios and ekpyrotic effects. This 
allows us to reproduce partially various models of modified gravity includ- 
ing EFG. If the Levi-Civita condition D|x=o=V is imposed on R we get 
modifications which are equivalent to that for f(R). In general, the cycling/ 
ekpyrotic models with nontrivial F are derived for certain nonholonomically 
induced torsion configuration. 



3.2 Models with effective anisotropic fluids 



To understand what is the difference between a usual f(R) theory and 
a Finsler type one with R(R) we have to consider physical equations with 
N-elongated partial derivatives and certain information of osculating ap- 
proximation ([2]). For simplicity, we shall study a model of effective locally 
anisotropic perfect fluid elaborated as follows. This allows us to write di- 
rectly generalizations of FLRW equations in N-adapted form, 



3H' Z = n 2 p, 3H' A + 2(e t H) 



-K 2 p, 



(22) 



where the energy-density and pressure of the locally anisotropic perfect fluid 
is defined in such a way that for R = R(R) 



P 



K- 2 {{d R R)- 1 

-K- 2 {(d R R)-' 



^R + 3He t (d R R) -3(e t H)}, 
^R + 2He t (d R R) + e t (e t (d R R)) 



(23) 



+ (e t H)}. 



We use such a N-adapted system of reference when @ and © are different 
in a form when pressure p of such a "dark" fluid is the same in all space- 
like directions (we remember that "tilde" is used for certain values encoding 
contributions from a nontrivial F). The equations of state (EoS) and EoS 
parameter in this Finsler model can be written/computed respectively 



P 



p + p and w = p/p, 



where 
i 



p : = - K - 2 {A(e t H) + e t (e t In \d R R\) + (e t In \d R R\) 2 } - He t In \d R R\} (24) 
has to be defined from a combination of FLRW equations, 

2e t H = k 2 x l p(H, e t H, e t (e t H), ...). (25) 
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It should be noted that in local coordinates such a system transform in 
a very cumbersome combination of functional and partial derivatives on t, 
with nontrivial N-coefhcients, which is quite difficult to be solved in ex- 
plicit form for some prescribed values Nf. We have to introduce additional 
frame/coordinate transforms and assumptions on H{t) in order to construct 
exact solutions for certain locally anisotropic "cosmic" functions (|24p . 

If the time variable is written as t(R(R)), it is possible to construct 
solutions of ([25]) following the approach developed in [25j 123]. For a new 
variable r(t,y(t)) = e t In \djiR(t)\, we transform the last equation into 

e t r + r 2 — Hr = 2e t H, for e t r = d^r — N^(t, y(t))d a r. 

To reconstruct the Hubble parameter for such locally anisotropic config- 
urations, we can consider typical examples with power-law or oscillating 
solutions. For instance, we show how we could generate oscillating Finsler 
configurations. We assume a particular behaviour when r = r + 1 r, for r(t) 
being the solution of 

d t r + r 2 -Hr = 2d t H (26) 

and a small value r is to be found from (neglecting terms ( r) 2 and 
N x ( V)) 

d t ( 1 r) + ( 1 r)(2r-Hr) = n (27) 

where n := N^d a (r — 2H). For n = and l r = 0, a solution of (|26p was 
found in [23j, when for r(t) = tq cosuiot it was reproduced 

f(R) = j dR exp[-r o; - 1 t(i?)]. 

In the case of Finsler modified gravity, the oscillating solutions of (|26p and 
(|27p can be expressed in the form r{t) = r$ cos wot + T\ cos uj\t. This which 
results in reconstructions of type 

R(R) = J dR exp { -roOJQ 1 sm[u) t(R)] - noj^ 1 sin[uJit(R)] } . (28) 

The function t(R) is an inversion of R(t) = 12H 2 + 6d t HE This corresponds 
to cyclic evolution reproduced both in f{R) and/or EFG gravity theories. 
The complexity of such solutions does not allow to obtain explicit forms of 
R(R) contributions. Nevertheless, we can distinguish from a "pure" f(R) 
theory and that with a mixed with F, when f(R,F) ~ R(R): in oscilla- 
tions of type ([28]) , there is the second term induced by off-diagonal terms 
summarized in n as a source of (|27p . 

8 A similar formula was derived in [26] for a model of osculating Finsler cosmology with 
weak anisotropy. 
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4 Tangent Lorentz Bundle Cosmology 



The goal of this section is to show how using R(R) and related f(R,F) 
cosmological scenarios on Lorentz manifolds V we can derive canonical 
models of Finsler gravity and cosmology on TV. Via Sasaki lifts of metrics 
([7]) we shall construct theories on total bundle spaces. 

4.1 Ekpyrotic and little rip Finsler cosmology 

Let us analyze scenarios of ekpyrotic / cyclic Universe \27\ [28] in our 
case derived from Finsler like modifications of GR. In this subsection, the 
constructions will be performed on /iTV components of tangent Lorentz 
bundles. Cyclic solutions solve the problems of standard cosmological model 
and provide a more complete theory. In ekpyrotic models, a scalar field is 
necessary for reproducing the cycling behavior and the f(R) gravity admits 
such phases in time evolution. Because the R(R) locally anisotropic gravity 
can be modelled as a f(R) theory, it is clear that both types of theories 
contain cyclic configurations. 

It is possible to reconstruct a canonical Finsler like model with a phan- 
tom phase and free of future singularity when it is generated a little rip 
cosmology similar to |29[ \30\ [3Tj . We can understand how such models may 
contain locally anisotropic modifications following arguments: Using the first 
formula in ()22|) . 3H 2 = K 2 p, the effective density function p (|23p and EoS 
parameter w = —1 — 2(e t H)/3H 2 , we find terms of type e t (djiR) and e t H 
with contributions from the N-connection coefficients. Big rip singularities 
(by definition) occur in a finite time t s when a(t) and the energy-density di- 
vergences. In f(R) models, this can be analyzed in local coordinate frames. 
Finsler modifications R(R) result in N-elongated partial derivatives ej and 
nonholonomic frames of reference. 

It is not clear if we can formulate a general condition for p and R(R) 
do avoid divergences. We can consider an enough but not sufficient con- 
dition that R(R) > but an additional investigation of contributions with 

for little rip. Using the reconstruction techniques from previous section, 
we shall prove that there are Finsler models of general acceleration, which 
do not contain future singularities and drive to a stronger growth in time 
than evolutions for big rip singularities. We can use the generating func- 
tion P(x) = exp[4/3e" x ], a = const, f3 = const, resulting in cycling factors 
of type (|19p . A similar function can be considered for different constants 
when we test the "cycling-accelerating" system for trivial N-coefficients, 
P(x) = exp[4/3e" x ], a = const, f3 = const. For small times t <C a, there are 
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reproduced de Sitter solutions and a Hubble parameter approximated to a 
constant. We get states with w < — 1, for a locally anisotropic phantom 
phase without big rip singularity. 

Let us show that a dissolution of bound structures (little rip) induced 
by a Finsler structure can occur. Using expansions of exponential functions 
into power series and solutions (1151) for certain values of R and/or R as a 
generating function for solutions of (|2ip . we compute three reconstruction 
models (see similar details in the beginning of Section V of [23] and references 
therein). The solutions depend on what type of scalar curvature we use, 
R(t) = 12H 2 + 6d t H, or R(t) = 12H 2 + 6e t H, and can be written in the 
form, 

f(R) = a 2 U + c 2y /4R/a* + 75) e V*A^+25/i6 

~ k x R + n 2 R 2 /a 2 + k 3j R 3 /" 4 + (29) 

R(R) = a 2 (ci + C*y/m/& + 75) 

~ kiR + k 2 R 2 /a 2 + k 3 R 3 /a 4 + 
f(R) = a 2 L + + 75^) eV«/i2o 2 +25/i6 

- kiR + k 2 R 2 /a 2 + K 3 ^ 3 /a 4 + 

where c\ = — 24exp (-39/12) and c 2 = 2\/3 are taken in order to obtain the 
same approximations if R — > R and /(...) — >GR; constants of type Ki, Ki, ... 
depend on c\,c\ and respectively on a,d etc. 

The above formulas with series decompositions distinguish three possi- 
ble cycling universes with rip evolution. In all cases, the quadratic terms 
on curvature allow to cure the singularities. But approximations are differ- 
ent: in the first case, we can obtain a recovering of GR form a "standard" 
f(R) gravity theory; in the second case, we model Finsler modifications as 
a /(...) theory without much information on Finsler generating function F; 
in the third case, the rip evolution starts from a Finsler modified spacetime. 
We suppose that such different cosmological scenarios can be verified ex- 
perimentally. For instance, for a Sun-Earth system with densities of type 
p = poe 2at ,p = poe 2at , and to = 13-7 Gyrs, we can obtain three different 
approximations for the time of little rip (decoupling) LR t, 

LR t = 13.7 Gyrs +29.9/a, LR t = 13.7 Gyrs +29.9/5, LR t = 13.7 Gyrs +29.9/a. 

For oscillating solutions of type ([28]) . we get possible resonant behaviour 
and shifting of decoupling etc. In general, models with f (R) , R(R) , f (R) 
possess different little rip properties. 
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4.2 Reproducing a canonical model of EFG 

A metric compatible Finsler model on TV can be completely defined by 
an action 

7( F R) 



S = Id x5 y\ \gijh, 



27? 



+ m L 



This would result in functional dependencies of type F R = F R(R) for 
a canonical scalar Finsler curvature ([9]). It is not clear how we could ex- 
tract information on a generalized gravitational constant (s) and matter field 
interactions via m L in extra "velocity /momentum" type dimensions. Nev- 
ertheless, we can encode such contributions into certain polarized cosmolog- 
ical constants derived for certain very general parameterizations of possible 
matter interactions with "velocity/accelertion" variables and their duals. A 
"recovering" of Finsler cosmology observational data on V to TV can be 
performed following such a procedure for two distinguished cases: 

1. Certain models of Finsler gravity can be associated to modified dis- 
persion relational 

u, 2 = c 2 [ % fc¥] 2 (i - l^ 2 .jJ l .J^/[^yW), 

where a corresponding frequency ui and wave vector fcj are computed 
locally when the local wave vectors ki — > pi ~ y a are related to vari- 
ables pi which are dual to "fiber" coordinates y a . These relations can be 
associated to a nonlinear quadratic element (see details in [32 | fl0 | [TT]). 

ds* = F 2 (x,y) « -{cAtf + < m {x k )y^[l + \ ^V^!^' ~ 1 + °(A 



r 



For physical applications related to "small" deformations of GR, we 
can consider that gij = (— 1, gj.(x k )) in the limit q — > correspond to 
a metric on a (pseudo) Riemannian manifold. 

2. Finsler variables with a generating function F(x, y) can be introduced 
in GR and f(R) modifications nonholonomic with 2+2 splitting. Such 
a F can be partially recovered in cosmological models using observa- 
tional data on /i-subspace. To understand possible physical implica- 
tions of theories on TV is important to construct exact solutions for lo- 
cally anisotropic black holes, brane trapping/warping and anisotropic 



9 for y l — dx' /dr, when cc*(r) is for a real parameter r; p%^ 2 i 2 (x) are parameterized 
by 3-d space coordinates with "hats" on indices running values i = 1, 2, 3 
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cosmological solutions \22\ 119] . The generic off-diagonal metrics, in 
both cases 1 and 2, are for Sasaki lifts ([7|). 



On tangent bundles, all fields depend on coordinates u a = (x l ,y a ) and 
the scalar curvature for the canonical d-connection D transform into R —> 
F R = R + R 0. Having recovered the spacetime metric hg = { gij(t) = 
gij(t, y(t))} for a cosmological model, we can construct up to frame transform 
a metric ([7|) , 



hg vg = gij(t) e* ® e j + h ab (x, y) e a ® e b 



(30) 



for certain coefficients h ab = e a a (x, y)e b (x, y)g a 'b'(t) and e a being deter- 
mined by a canonical N-connection N?(x,y) induced via a chosen F(x,y), 
see footnote EJ For cosmological models, the coefficients of v -metric can 
be transformed via frame transform to h a b(t,y). In general, such mod- 
els are inhomogeneous on fibre coordinates y a . Extending the scalar fields 
x(x l ) — > x{ xi iV a )i we can construct values of type (fTT|) . F P[x( F R)] > 
and F Q[x( F R)}, when F R(R) = F P[ X (R)]R+ F Q[ X (R)}- 
Actions of type (|20p can be generalized in the form 



S 



d 4 x5 4 y^\g ij \\h ab \ 



-L F R(R, X )+ ™ L 



where 5y a = e a ([5]). We do not consider a factor e~ 2x in this action because 
we shall work with another type of conformal transforms when the action is 
re-defined in Finsler generalized Einstein frames in order to remove strong 
coupling. On TV, a quintessence locally anisotropic action is postulated in 
the form 



S 



/' 



d xd yy\gij\\h ab \ 



1 



R--u(x)(D a x)(D a x)-U( X ) 



where g Q/3 = F P( X )g a f3, w = 12 ^^Pj / F P,U= F Q/( F P) and D 
is the conformal transform of D. For cosmological models, there are used 
(in our case, Finsler analogs) the Jordan frames with 



a(t) 



7 P( X (t))a(t), for dt = J F P( X (t))dt, 



and locally anisotropic configurations x(^y a )i where the fiber coordinates 
y a can be parameterized y a (t). In EFG, the evolutions with a(t) and a(t) 
are with respect to N-adapted frames (jl]) and ((SJ). 
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We reproduce cycling universes with rip evolution of type (j29|) prescrib- 
ing^| 

u(x) = 43/3 2 e 2Sx and U( X ) = -6S 2 (3 + 4/3e s *)(3 + 8f3e Sx ) exp(-4/3e Sx ). 

Keeping only terms with i-evolution and e at ~ 1 + at + 0(t 2 ) for small t, 
the solution for scaling factors are respectively constructed 

^ _ ^ /"2^exp[Sx] 

a(i) = aoe[6(/3e ax + at)], for ao = const and i = / e 2 /z; (31) 

a(i) = ao? 6 exp[6a/3e _2/3 i], for ao = const and = S/3e _2/3 + 6/t. 

Exact solutions in EFG with such scaling factors will be constructed in next 
section. Here we observe that the Universe with Einstein-Finsler frames 
describe both a type of (initial) big bang singularity and a super-accelerating 
evolution. This is a manifestation of scenarios with little rip determined by 
possible locally anisotropic character of gravitational interactions on TV. 
It seems that singularities can be removed in Jordan frames adapted to N- 
connections. We conclude that Finsler configurations extended on tangent 
Lorentz bundles may result in dissolution of bound structures of certain 
classes of FLRW models originally defined in GR and then extended to 
EFG. 

4.3 Cosmological solutions in Finsler gravity 

On TV of a Lorentz manifold V, the Einstein-Finsler equations (jlOp for 
the canonical d-connection D are for a 8-d spacetime endowed with non- 
trivial N-connection structure. Such systems of nonlinear partial derivative 
equations (PDE) can be solved in very general forms using the anholonomic 
frame method [22\ [32l [2Tj . Locally anisotropic and Finsler like solutions in 
4-d models of gravity were constructed in Refs. (Ill 1191 \12\ . In this subsec- 
tion, we provide several examples of 8-d exact solutions which possess cyclic 
/ ekpyrotic and little rip properties. 

4.3.1 Decoupling of EFG cosmological equations 

We label local coordinates x % = (x 1 = r, x 2 = 9),y ai = (y 3 = t, y 4 = ip); 
y a = (y a2 ,y a3 ), where indices run respective values a\ = 3, 4; 02 = 5, 6; 03 = 

10 we extend on TV similar formulas (59) in [23]; the coefficients are redefined with 
"hats" and fixed in such a form on ft-subspaces we get little rips as in f(R); nevertheless, 
the parametric dependence is modified on velocity like coordinates and parameters for 
EFG 
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7,8. We consider (x l ,y ai ) as coordinates on a 4-d Lorentz manifold V, and 
the coordinates y a as fiber coordinates in TV. This reflects a conventional 
2 + 2 + 2 + 2 splitting of coordinates which will give us the possibility to inte- 
grate the gravitational field equations "shell by shell" increasing dimensions 
by 2. The FLRW cosmological solution can be written in the form 



hg=a 2 (t) I - — — + r 2 d9®d6j - dt®dt + a 2 (t)r 2 sin 2 6d(p®d(p, (32) 

with a = 0, ±1. This metric is an exact solution of the Einstein equations 
with a perfect fluid stress-energy tensor, T a a = diag[—p, —p, p, —p], where p 
and p are the proper energy density and pressure in the fluid rest frame. For 
simplicity, we consider as a "prime" a /i-metric ansatz containing conformal 
transform of (132f) . multiplying on a~ 2 , with a = 0, and in Cartesian coordi- 
nates, hg = g i dx t ®dx' 1 +h a dy a 0dy a , where g~i = 1, h^(y 3 ) = — a~ 2 (t), = 1. 
To construct 8-d cosmological Finsler like solutions, we shall use the ansatz 

g = r] t {x l )gdx l ®dx l + i lai (x\t)h ai e ai ®e ai + 

h a2 {x k , t, y 5 )e a ^e a2 + h a3 (x k 7 t, y a * , y 7 )e a3 ®e a \ 
e ai = dy ai +N? 1 (x k ,t)dx\ (33) 
e" 2 = dy a2 +N^(x k ,t,y 5 )dx l +N^(x k ,t,y 5 )dy a \ 
e a3 = dy aa +N? 3 (x k ,t,y a2 ,y 7 )dx l + 

N% (x k , t, y a - , y 7 )dy^ + N a a * (x k , t, y a * , y 7 )dy a - , 

where h ai = rj ai (t)h ai are defined by polarization functions r] ai = 1 + 
e Xa± (x l ,y 3 ), e < 1. This class of generic off-diagonal metrics depend on time 
like coordinate, t, and on fiber-like ones, y a2 and y 7 , and possess a Killing 
symmetry on d/dy 7 because the coefficients do not contain the coordinate 
y 7 . In Cartesian coordinates, the fe-part of ansatz ([33]) is modified by rfc = 
1 + exi(x l ), when g t = rjigi = e^^^O 

We do not have explicit observational/experimental evidences to deter- 
mine what kind of sources Y^s may have physical importance for models 
of matter field interactions in the total space TV. From a formal point of 
view, we can extend a geometric/ variational formalism for deriving energy- 
momentum tensors on V (for instance, for scalar, spinor, gauge etc fields) 
to construct similar values using Sasaki lifts for metrics and adapting the 
constructions to N-elongated frames ([!]) and ([5|). For simplicity, we shall 



11 We note that it is possible to construct more general classes of solutions without Killing 
symmetries, nonhomogeneous cosmological metrics etc. Such an ansatz is a natural one 
when on the 4-d base spacetimes the cosmological metrics depend only on t and can be 
modified by 3-d space coordinates and/or several fiber type coordinates. 
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approximate such possible "extra velocity" contributions by matter interac- 
tions with an effective cosmological constant A when T% = A5p. Then we 
shall see what kind of 8-d generic off-diagonal cosmological solutions would 
possess cyclic, ekpyrotic and/or little rip properties in the h-part of the total 
metric. As a matter of principle, such an effective cosmological constant can 
be "polarized" and depends on base and fiber coordinates. Nevertheless, it 
can be transformed into a constant value using re-definition of generating 
functions for various classes of solutions as we shall prove below. There 
are also possible certain mechanisms of Finsler brane trapping/warping of 
extra "velocity" type coordinates when N-adapted geometric constructions 
are considered [22j [32] . 

The gravitational field equations (fTUj) for the ansatz (f33j) decouple in this 
form: 

In 4-d, we get a system of nonlinear PDEs which with respect to INT- 
adapted frames is written 

dt^ + d 2 ^ = A, (34) 

</>*(lnN)* = A/i 3 , (35) 
PNf + Oi = 0, (36) 
(iV/)** + 7(A?)* = 0, (37) 

with = d t (f>,di = 83.1,8^ = 3\ !, where the coefficients 

7 = (In |/i 4 | 3 / 2 - In \h 3 \)*, cti = h\di<t> and = hlf, (38) 

are determined by /13 and /14 via generating function 

<f>(x i ,t) = hi\2(]ny/\h^)*\-hiy/\f^\. (39) 

The equations on the first 2-d shell (with coordinates y 5 and y 6 ) are 

(8 5 V)d 5 (m|/t 6 |) = A^, (40) 
1 0N? + l ai = 0, + V = 0, (41) 

dUNf)+ Vs(A?) = 0,5 5 2 5 (iV a 6 1 )+ V 5 «) = 0, (42) 

for coefficients 

x 7 = fc(ln|/i6| 3 / 2 -]iiN), 1 a l = (8 5 h 6 )(8 l 1 4>), (43) 
l a ai = (dsheXda, 1 ®, 1 /3 = (d 5 h 6 )(d 5 V), 

determined by /15 and h,Q via 

V(x fc ,t,y 5 ) = ln|2«9 5 (ln^i)| -lnyW- (44) 
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In similar form, using coefficients h 7 and hs (and third "shell coordinates" 
y and y 6 ) we obtain 

(&7 2 cj>)d 7 (\n\h 8 \) = Ah 7 , (45) 

2 $Nj + 2 a t = 0, 2 pN 7 ai + 2 a ai = 0, 2 /3iV a 7 2 + 2 a a2 = 0, (46) 

d 2 7 {Nf)+ 2 1 d,{Nt) = 0,d 2 7 (N 8 ai )+ V 7 «)=0, (47) 

5 7 2 7 «)+ 2 79 7 «) = 0, 

for coefficients 

2 7 = £>7(lnN 3/2 -lriN), 2 a i = (d 7 h 8 )(d i 2 <P), (48) 
2 a ai = {d 7 hs)(d ax 2 4>), 2 a a2 = (drhs^da, 2 ^, 2 (3 = (d 7 h 8 )(d 7 2 </>), 

generated by 

2 (j ) (x k ,t,y a \y 7 ) = ln|2d 7 (ln^7^)| - In yW. (49) 



We can see that the system of equations ([35]) - ([57]) is similar, respectively, 
to (130]) - (132]) and ([35]) - (137]) , Such equations can be integrated consequently 
by adding additional dependencies on next shell coordinates. 



4.3.2 Generating off diagonal cosmological Finsler solutions 

The fr— metric is given by e^^ x ^dx i <g> dx\ where tp(x k ) is a solution of 
([33]) considered as a 2-d Laplace equation ([33]) . It depends on effective 
cosmological constant A. 

We can integrate the system ([35]) and ([39]) . for <fi* ^ 0. Such a condition 
can be satisfied by choosing a corresponding system of frames/coordinates; 
it is possible to construct solutions choosing <p with (p* = 0, as particular 
cases (for simplicity, we omit such considerations in this paper). Defining 
A := (In \ h4\y and B = y/\h%\, we re-write 

4>*A = KB 2 , Be 1 ^ = 2A. (50) 

If B 7^ 0, we get B = (e^)*/2A as a solution of a system of quadratic 
algebraic equations. This formula can be integrated on dt which results in 

h 3 = %(1 + (^)*/2Av^^) 2 . (51) 

Introducing this /13 in (|50p and integrating on t, we get the coefficients 

h A = O/Hexpt^-H, (52) 
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for an integration function /i4 = /i4(x fc ). We can fix °h a = h a as in (|33[) . 

Having defined h a we can compute the N-connection coefficients as so- 
lutions of ([36]) and (integrating two times on t) ([37]) . 

Wi = -d i( j)/4>\ (53) 
nk = l nk + 2 n k J dt h 3 / (y^\fu\) 3 , 

for integration functions 1 rik{x l ) and 2 rik(x l ). 

Introducing solutions (|5ip , (|52p and (|53|) for the /i-metric of ansatz (|33p , 
we get an quadratic element for nonhomogeneous 4-d cosmologies, 

ds 2 = e^(dx') 2 + h 3 (l + -^-£=? [dt - ^dx f + (54) 

2A V /|A 3 | <?* 

hexp[^W + 2 n k J dtj^L=^)dx 1 ] 2 . 

The solutions depend on generating functions (f>(x l ,t) and tp(x k ) and on inte- 
gration functions l nk(x k ), 2 rik(x k ). We approach the FLRW metric (|32p if 
we chose such values that r\i = l+exi{x % ) — > 1 and n ai = l+ex ai {x l ,y 3 ) — > 1 
and the N-connection coefficients vanisrl^l. This class of modified Finsler 
like spacetimes (effectively modelled for 2+2 decompositions) are character- 
ized by a nontrivial torsion field determined only the coefficients of metric 
(and respective N-connection). So, extra-dimensional Finsler contributions 
can be via off-diagonal extension of solutions, nonlinear polarizations of 
physical constants and metric coefficients and induced torsion. Such metrics 
(being generic off-diagonal) can not be diagonalized via coordinate trans- 
forms because the anholonomy coefficients in are not zero for arbitrary 
generating/ functions. 

Metric of type (|54p can be constrained additionally in order to construct 
exact solutions for the Levi-Civita connection V. We have to consider solu- 
tions with 2 nk = 0, di( l nk) = dk{ 1 rii) and Wi = —di^/cj)* and subjected 
to 

w* = &i In \h±\,diWj = djWi,n* = 0, 

see details how to solve such equations in [201 1221 ESJ [21] (this is possible, 
for instance, for separation of variables). Even for solutions with V, we 

12 Here we note that the FLRW solution is for perfect fluid stress-energy tensor and 
not for a cosmological constant. Nevertheless, using generating and integration functions 
for coupled nonlinear systems we can define limits to necessary type diagonal metrics via 
corresponding nonholonomic constraints and frame transforms. 
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get only effective Einstein spaces with locally anisotropic polarizations and 
off-diagonal terms induced from Finsler gravity. 

Let us consider extra shell fiber contributions. The equations (|40p - ()42p 
and ([4"5]) - (|4"T|) can be integrated following above procedure but for corre- 
sponding coefficients, ([4~3]) and ([Ml), and generating functions, (|4"4"]l and (fH)]) . 
Such a class of 8-d solutions are parameterized by the quadratic element 



ds 2 = e^(dx*) 2 + h 3 (l + dt ^ e p_ f [dt - ^dx 1 } 2 

2A V /| h 3 \ 

h 3 



(55) 



dt- 



+ °ft 5 (l + 



K exp[ 



d 5 e 



f [dy 5 - ^Mdo 



o2 1 <t> 



8A 



-][dy 6 + ( 1 iV fe 6 + 2 7V£ / dy 



)dx l 



(VM) 3 

ftg 



oil2 



(vW) 3 



+( ^ + ^< / dy 



)dy 



a 1 1 2 



°ft?(l 



)cte l 



3 /i8 exp[ 



3 7 ( 2 0) 



a 2 l2 



.2 ^ 



8A 



-][d J/ 8 + ( 1 iV|+ 2 AT| / dy 7 



)dx l 



+( x < + 2 < / dy 7 



ft? 



)dy a > + ( X + 2 < 



ft? 



)dy a \ 



where there are used the integration functions 1 A^, 2 N^ depending on 
(x\y 3 = t); °h ai , 1 iV«, 2 iV fe 6 , ^ 2 JV* depending on (x\y a \y 5 ); 
°h a2 , x Nl, 2 N 8 k , l N^, 2 N^, 1 N^, 2 iV a 8 2 depending on O^^y 7 ). 
The generating functions are 

We can fix such configurations which model certain scenarios in cosmology 
and/or gravitational and matter field interactions. 



4.3.3 Extracting realistic cosmological configurations in EFG 

Metrics of type ()55[) define off-diagonal exact solutions for 8-d locally 
anistoropic generalizations of FLRW cosmology. In general, it is not clear 
what kind of physical significance such nonhomogeneous solutions may have. 
Playing with values and parameters in generating and integration functions 
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and source A, we can mimic different scenarios and compare them, for in- 
stance, with modifications derived for other types of gravity theories. 

We note that we can generate off-diagonal solutions depending only 
on t and fiber coordinates (y a2 ,y 7 ) if via frame/coordinate transforms and 
a corresponding fixing of generating/integration functions when the coeffi- 
cients of ()55[) do not depend on coordinates (x l ). In such cases, a series of 
coefficients N® = and, the integration functions are of type (p(t),dt4> 7^ 
0; V(t,y 5 ),c>5 l <j) 7^ 0; 2 (j)(t,y 5 ,y G ,y 7 ), dj 2 <j) 7^ 0. The nonholonomic/ non- 
linear gravitational dynamics on fiber variables may mimic scalar and other 
matter field interactions in observable 4-d spacetime. It is possible to apply 
various trapping/warping mechanisms like in brane gravity [221 E2] or using 
osculating approximations in Finsler gravity. 

Introducing a cycling factor a(t) = exp[Hot + b(t)] (fl"9l) instead of a(t) 
in /i3(y 3 ) = —a~ 2 (t), we can model such a cycling scenarios with respect 
to N-elongated frames (j4]) and © when the nontrivial coefficients AT" are 
functions on some variables from the set (t, y 5 ,y e , y 7 ). For small values of N 
(assuming that fiber gravity results only in small corrections), we can put 
certain symmetries and boundary conditions on a subclass of solutions (|55p . 
which possess cycling behavior and may limit in 4-d the FLRW solution 
(|32p but contains also certain nontrivial anisotropic polarizations caused 
by possible Finsler like interactions. For anisotropic spacetimes, we have to 
recover not a f(R) theory and certain "exotic" states of matter but choosing 
a corresponding set generating/integration functions. 

To reproduce little rip evolution by off-diagonal Finsler metrics we can 
introduce a scaling factor a(t) = ooe[6(/3e ax + at)] (f3Tj) instead of a(t) in 
hs(y 3 ) = —a~ 2 (t), or try to get such a term in /13 and /14, for (|55p . All 
geometric constructions will evolve with respect to certain nonholonomic 
frames of reference, with small corrections by N-coefficients. 

Finally we note that a solution of type (|55f) is not written in well known 
Finsler form (|7|) with standard Hessian g a f, ([1]) and N-connection structure 
iV". Nevertheless, we can related the coefficients of both representations 
using frame transforms of type g a iy = e a a ,e b b ,g a i } if a fundamental Finsler 
function F is taken from certain experimental data or fixed following certain 
theoretical arguments. It is not convenient to find exact solutions working 
directly with data yg^jN^j because the gravitational field equations con- 
tain in such cases forth order derivatives on F. To derive a general geometric 
method of constructing exact solutions is possible for ansatz of type (|55p with 
general dependence on certain classes of generating/integration functions. 
Then constraining the integral varieties of solutions and after corresponding 
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frame transforms and distorting of connections ([8]) we can derive exact solu- 
tions for necessary types of connections, symmetries, cosmological evolution 
behavior which may explain cyclic and little rip properties etc. 

5 Discussion and Conclusions 

The procedure for explicit reconstruction of modified gravity using cos- 
mological observation data [231 [24"1 [25] was generalized for a comparative 
study of f(R) and Einstein-Finsler gravity (EFG) theories. We have shown 
how such models can be extended around general relativity (GR) action 
and various classes of cosmological solutions. The conclusion on existence 
of cyclic evolution scenarios driven by Finsler fundamental functions was 
made using locally anisotropic variants of scalar-tensor gravity theories. 
Nonholonomic constraints and off-diagonal metrics may mimic matter field 
interactions and evolution processes. Geometrical actions, off-diagonal non- 
linear dynamics and the conditions usually considered for inflation models 
were shown that may lead to cyclic universes and ekpyrotic effects. 

We proved that both f(R) and EFG theories encode scenarios with lit- 
tle rip universe with dark energy represented by nonholonomically induced 
non-singular phantom cosmology. So, such classed of theories are consistent 
with observational/ experimental data and present a realistic alternative to 
the ACDM model. Various types of non-singular super-accelerating and/or 
locally accelerated universes can be reconstructed in modified gravity the- 
ories. The solutions presented here can be extended to Finsler like brane 
and/or Hofava-Lifshitz theories [22^132] . Such theories modelled as GR plus 
corrections pass various local tests. 

It is a very complicate technical problem to construct generic off-diagonal 
solutions in GR and modifications. We generalized our anholonomic frame 
method |20[ \2l\ of constructing exact solutions in a form which would allow 
to decouple and solve in off-diagonal forms the gravitational field equations 
for modified theories. Following such an approach, it is possible to gen- 
erate anisotropic off-diagonal cosmological scenarios with dependence on 
"velocity- momentum" type coordinates y a , see metrics (155p . Such models 
can be considered as low energy limits of Finsler stochastic metrics of type 
(3.4), (3.10) and (4.18) in [16], when y a can be associated to momentum 
transfers by quantum-gravitational fluctuations in the space-time and D- 
brane/-particle foam. It is possible to provide a microscopic background 
following this approach. Our class of solutions allow us to model cosmologi- 
cal evolution scenarios with acceleration and dark energly/matter effects in a 
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form adapted to the nonlinear connection structure (with is very important 
in Finsler gravity) . For various types of Finsler gravity and cosmology theo- 
ries with y a treated as extra dimension coordinates, the physical meaning of 
such variables is model dependent. For instance, we can not perform com- 
pactification of velocity type coordinates because there is a constant speed 
of light and nontrivial N-connection structure. It is important to find solu- 
tions with trapping/warping effects. Having constructed certain physically 
important Finsler gravity /cosmology models we can perform the osculating 
approximation y{x) to metrics of type ([2]). This allow us to compare obser- 
vational data for different modified theories of gravity and GR which in all 
case are considered in real/effective 4-d spacetimel 13 ! 

Finally, we not that in this paper the reconstruction procedure is involved 
in a new form when the cosmological models are determined by certain 
classes of generating/integration functions. The cyclic/ekpyrotic/ little rip 
universe scenarios are possible in EFG but it remains to understand how 
natural and realistic are solutions with locally anisotropic effects on tangent 
bundles on Lorentz manifolds. This is one of the scopes of our future works. 
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13 There were elaborated various geometric and physical models of Finsler spacetimes 
and applications Finsler methods in modern gravity and cosmology, see |13l 1141 1151 1161 1171 
18, 19 20 21 22, 33 34 and references therein. A series of works are with metric noncom- 
patible Finsler type connections, or without explicit assumptions and physical motivations 
on Finsler nonlinear and linear connection structures. The main physical problems of such 
works with nonmetricity fields are that we are not able to define analogs of Finsler type 
spinors, Finsler-Dirac equations and well-defined "standard" conservation laws. There 
are not "simple" analytic methods for constructing physically important exact solutions 
for Finsler type black holes, brane configurations and to formulate a recovering proce- 
dure from cosmological data became very problematic etc, see critical remarks in Refs. 
[9j 1111 1101 [8j [12]. Perhaps, most closed to general relativity and standard physics are 
the (generalized/modified) Finsler like models constructed on tangent bundles to Lorentz 
manifolds when there are defined natural lifts of geometric/physical objects from Einstein 
spaces to Einstein-Finsler analogs. In such cases, there are canonical metric compati- 
ble Finsler connections adapted to the nonlinear connection structure and an osculating 
approximation can be performed. Following such an approach, it was formulated a self- 
consistent axiomatic formalism [11] which is very similar to that for the general relativity 
theory. It is possible to elaborate certain renormalizable models of Finsler like quantum 
gravity [32] and recovering procedures from cosmological data, to find exact solutions 
[22] [201 [21] etc. 
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